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ABSTRACT 

Light bending due to strong gravity has recently been invoked to explain variability and flux 
correlations between different bands in some accreting black holes. A characteristic feature 
of light bending is reflection-dominated spectra, especially if photon sources lie in the deep- 
est parts of the gravitational potential within a few gravitational radii of the event horizon. 
We use the spectrum of the hard X-ray background in order to constrain the prevalence of 
such reflection-dominated sources. We first emphasize the need for reflection and explore the 
broad-band properties of realistic spectra that incorporate light bending. We then use these 
spectra, in conjunction with the observed 2-10 keV AGN distribution, evolutionary and ob- 
scuration functions in order to predict the hard X-ray background spectrum over 3 — 100 keV, 
and provide limits on the fraction of reflection-dominated objects, dependent on the height of 
the photon sources. Our results allow for a cosmologically-significant fraction of sources that 
incorporate strong light bending. The luminosity function based on intrinsic flare luminosities 
is derived and implications discussed. We discuss prospects for future hard X-ray missions 
such as NeXT and Simbol-Xthat can image such sources as well as confirm the precise spec- 
tral shape of the background near its peak, important for constraining the universal relevance 
of light bending. 

Key words: X-rays: galaxies - X-rays: extragalactic background - galaxies: active 



1 INTRODUCTION 

Refle ction and scattering of X-ray photons on optic ally-thick media 
(e.g. [Pounds et all 1 19901 : lOeorge & FabianI 1 1 99 ih are believed to 
produce features such as the Fe fluorescence emission hnes and the 
broad 'Compton hump' in the spectra of accreting Galactic black 
holes (BHs) as well as active galactic nuclei (AGN), allowing the 
geometry and rotational dy namics of matter accreting onto the cen - 
tral BH to be probed (e.g. Miller et al."2004'; Fabian et al.'2002a|). 
[Mini utti & Fabian (2004, see also Fabian & Vaughan 2003) have 
recently shown that incorporation of general relativistic light bend- 
ing effects on photon propagation can have dramatic consequences 
for the spectra of accreting sources with peaked central emissivi- 
ties. For isotropic photon sources (such as magnetic reconnection 
events/flares, a hot electron corona or the base of a jet) that lie 
within a few gravitational radii (_Rg = GA/bh/c^) of the BH, the 
cross-section of photon geodesies intercepting the accretion disk 
(or lost down the event horizon) can be very large, as compared to 
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sources at much higher heights (h). For constant luminosity photon 
sources, the primary flux observed at infinity can vary dramatically 
with changing source height, while the reflected flux component re- 
mains relatively unchanged, resulting in reflection-dominated spec- 
tra for low h. This model has been shown to be capable of ex- 
plaining the apparently dis-jointed variability patterns and flux-flux 
correlations observed between different bands in some Seyferts, 
Narrow Line Seyfert i s, as well as Galactic BH candidates (e.g . 
Miniutti et al.ll2007allbl l2004l : iPonti et al.ll200^ : iFabian et al.ll200g. 
2004lRossi et al.ll2005l) . 



The general importance of the light bending model is currently 
unclear The gravitationally blurred emission lines and peaked cen- 
tral emissivities of the sources for which light bending has been 
proposed suggests that the emission is being generated, or repro- 
cessed by the disk, within a few _Rg of the BH. If the disk is as- 
sumed to extend down to the last stable circular orbit (and assum- 
ing no significant emission from radii smaller than this orbit), then 
radii of less than 6Rg (the Schwarzschild limit) are inferred. Since 
smaller stable orbits are possible in the rotating Kerr metric, these 
BHs are inferred to be spinning. Loose constraints on the spin pa- 
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rameter can be placed by studying the amount of gravitational red- 
shift implied by the Fe line, and suggest that the BH spin is prob- 
ably rapid. It is also in the inner parts of the gravitational potential 
well that light bending effects are strong. Other increasing evidence 
also indicates tha t most super-mass i ve black holes (SBHs) must be 
rapidly spinning jWang et al.ll2006l : lElvis et aljbood ; iFabian et af] 
l2002lJ : l\^onteri et alj|2005h . If true, then light bending should be 



considered seriously in a cosmological context. Future X-ray mis- 
sions such as XEUS, NeXT and SimboJ-Xwith large effective areas 
for imaging and spectroscopy over wide energy bands will search 
for the characteristics of light bending such as strong, blurred emis- 
sion lines and high reflection fractions. On the other hand, it should 
also be possible to place some constraints on the fraction of sources 
with reflection-dominated spectra by using the hard X-ray back- 
ground (XRB) radiation spectrum, since this background is simply 
the integrated emission of AGN spread out in redshift. This is what 
we aim to do in this paper 

The hard slope of the X-ray background (XRB) radiation and 
its pe ak around 30 keV have been well-known for several decades 
now l lKinzeretal .1119781 : [Marshall et all 1 19801 : iFabian & Barcon j 
I1992I and references therein). Recent surveys with ^and XMM- 
Newton have proven that the 2-10 keV XRB (and probably that up 
to ~ 100 keV) is dominated by the integrated emission from a com- 
bin ation of obscured and uno bscured AGN (see, e.g., the review 
by iBrandt & Hasii^ l2005l. and references therein). XRB syn- 
thesis models have had much success in simultaneously satisfying 
several observational constraints including the slope and normal- 
ization of the XRB spectrum, the number counts of AGN resolved 
below 10 keV, as well as the A GN redshift distribution observed 
in deep field o bservations (e.g., IComastri et alj[l99^ : iGilli et all 
I2OOII : lGandhr& Fabian 2003). Recent models are also able to fit 
the observed number counts in the mid-infrared and the sub-mm, 
as well as the fracti on of obscured sources as functions of luminos- 
ity and/or redshift jTreister & Ugvll2005l : iBallantvne et al.l I2OO6I : 



iGranato et all2006l : 



Gilli et al.ll200l) 



The characteristic hard slope and 40 keV break of the XRB 
suggests the presence of at least some reflection in AGN spectral 
energy distributions (SEDs), and the above models typically use a 
reflection fraction {R) ~ 1. Yet, few firm constraints on the high en- 
ergy spectral shape of AGN exist, in comparison to those below 10 
keV. While studies ^i\hBeppoSAX, and now with Swift and INTE- 
GRAL, have begun to probe large samples of AGN over the energy 
range corresponding to the peak of the XRB, these are still sensi- 
tive mainly to the bright and loc al population, and resolve only ~ 1 
per cent of the relevant XRB ( Dadina 2007: Ma rkwardt et alj2005l : 
iBeckmaim et al. 2006). A handful of sources with high R have 
recently been confirm ed with the S uzaku satellite (MCG -6-30-15 
iMiniuttietalJl2007all : NGC 4992 jComastri et al.ll2007ll ). but the 
abundance of such sources is ill-constrained, though th e existence 
of a much larger population is hinted at by the data dUeda et al] 
l2007h . 

In this paper, we use realistic reflection-dominated SEDs in 
conjunction with the observed 2-10 keV number counts and evolu- 
tion of AGN in order to constrain the fraction of sources in which 
light bending can be an important effect. We do this by fitting to the 
flux and spectral shape of the hard XRB spectrum as measured by 
the HEAO satellite, allowing for moderate adjustment of the AGN 
number density (and, in particular, the fraction of Compton-thick 
AGN) within the current XRB normalization uncertainti es. As our 
base SEDs, we use the ones of lSuebsuwong et al] ( l2006l hereafter 
S06), who have carried out detailed Monte Carlo investigations that 
incorporate light bending to quantify changes in observed spec- 



tral properties, with the variables being inclination angle, source 
height, source size (either point-like or ring-like) and source dis- 
tance from the rotation axis, assuming the proper Kerr space-time 
metric and no emission from the plunging region between the in- 
nermost marginally stable orbit of the accretion flow and the event 
horizon. Th e specific ang ular momentum of the BH is taken to be 
a = 0.998 ( iThomelf 19740 . implying that the innermost stable orbit 
lies at a distance of 1.23 Rg from the BH. For simplicity, we con- 
sider only the case with on-axis point-like flares at height h above 
the accretion disk. Finally, these base SEDs are obscured in order 
to match the column density obscuration function determined from 
surveys below 10 keV. 

We note that our work places rough constraints on the frac- 
tion of reflection-dominated sources in general, irrespective of 
the physical mechan isms that generate their spectra. For instance, 
iMerloni et all ( I2OO6I) discuss a model wherein radiatively-efficient 
sources are subject to highly inhomogeneous accretion flows due 
to radiation-pressure related clumping instabilities. They show that 
reflection-dominated spectra can result, provided that energy re- 
lease is highly localized and the sky covering factor of the cold 
phase of the inhomogeneous flow is large enough. Their model 
can reproduce the dis-jointed variability of the primary power-law 
and the reflected component, in agreement with the light bending 
model. Distinguishing between these various models is not the pur- 
pose of this paper. 

We assume Ho = 70 km s"^ Mpc"\ Qm = 0.3 and A = 
0.7, where required. 



2 THE IMPORTANT OF REFLECTION FOR THE XRB 

Before detailing the effect of light bending, we consider the ne- 
cessity of reflection in order to fit the XRB spectrum. The pos- 
sibility that reflecti on could exp l ain th e XRB spectral hump was 
first pointed out bv [ Fabian et alj ( 1 19901 . and commented upon by 
IZdziarskiet"aLlll993l) r and reflection is now typically included in 
XRB synthesis models, as mentioned in the previous section. There 
is renewed interest in this issue, however, due to recent results from 
sensitive, hard X-ray detectors (especially SwiftfBAT; see below) 
that are beginning to place stricter observational constraints on this 
parameter. We thus provide some detailed discussion on the issue, 
by asking whether it is possible to make do without including re- 
flection. 

Fig [B (top) shows a XRB model spectrum assuming the 
generally-accepted 2-10 keV X-ray lu minosity function (XLF) 
and obscuring column (A^h) function of ' Ueda et al ]|200l here- 
after U03), including a fraction of Compton-thick (logA^H=24- 
25) AGN determined from the extrapolation of the A^h func- 
tion beyond logA ^ H=24. The broad-band SEDs used are those of 
iGandhi & Fabiaril tOOl hereafter GF03) with parameters very 
similar to those of U03, but without any refl ection co mponent in- 
cluded. The observed XRB spectrum is that o f iGruber et al. ( 199^, 
as scaled by U03. More details on the XLF, the SEDs and other as- 
sumptions (most are used again in this paper) will be given in the 
following sections. Since we only aim to emphasize the effect of 
reflection in this section, minimal details are stated. The total mod- 
eled background (thick solid line) underestimates the observed one 
by ^ 15, 25, 39 and 37 per cent at 5, 10, 30 and 50 keV respec- 
tively. Without reflection, the spectrum is not 'peaky' enough (see 
also Fig. 18 of U03). 
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2.1 An additional population of peaky sources? 

A natural first question to ask is wliether some additional popula- 
tion of sources (eittier known or unknown) could compensate tiie 
shortfall in the background. This would have to be a population 
with peaky (reflection-free) broad-band spectra in order to recover 
the deficit seen in the figure. The complete absorption of radiation 
below 10 keV in the rest-frame for Compton-thick AGN also results 
in net spectra with narrow and peaky shapes. While recent 2-10 
keV have constrained the number densities of Compton-thin AGN 
well, these very surveys probably missed many Compton-thick 
AGN. So how large an increase (if any) in the number of Compton- 
thick sources would suffice to recover the shortfall in the XRB? 
Fig. [T] shows (dashed lines) that the spectrum produced by a uni- 
form 4-fold increase in the Compton-thick population (logA'H=24- 
25) at all redshifts is required to recover the deficit at all energies. 

There is, however, no compelling evidence to support the ex- 
istence of such a large populati on, either in the optical or in hard X- 
rays (at least at low redshift : iRisaliti et alJ ll993 : lMarkwardt et al] 
l2005l : lBeckmann et alj2006l) . Additionally, the total BH mass den- 
sity with this increased population of sources is 7.8-19 xlO^ 
M0Mpc^'^, assuming that BHs grow by accretion with a radia- 
tive efficiency 77=0.! and a Bolometric:2-10 keV luminosity cor- 
rection factor with a range of 20-50 (Vasudevan & Fabian 2007, 
submitted). This exceeds t he upper limits sugg ested by other re- 
cent studies ( lFabianll2004l :l Marconi et al. 1 12004 "). No other simple 
modification of the model (e.g. an increase in the Compton-thick 
fraction with redshift) produces an acceptable solution; we do not 
investigate more exotic evolutionary scenarios. 



2.2 Including spectra with harder photon indices? 

lGimetal.l ( l2007l) have shown that the modeled XRB spectrum can 
be hardened, especially near its 30 keV peak, by accounting for 
a dispersion of power-law photon indices. Could this be another 
way to compensate for the hard X-ray flux that is missed when 
no reflection is included? Fig. [T] {bottom) shows such a fit to the 
XRB. Choosing a realistic gaussian dispersion or ~ 0.2 (for an 
average photon-index To = 1.9) increases the 30 keV flux by 16 
per cent (dotted line) above the assumption of ar = (solid line), 
but still underpredicts the observed XRB. A much larger dispersion 
of or ~ 0.5 (dashed line) is required to harden the slope enough 
to match that of the XRB below 10 keV. This then results in a large 
excess of the background flux above 10 keV, eve n for high energy 
exponential cut-off (Scut) of 200 keV chosen bv lGilh et al] llool 
already lower than values of 360 and 500 keV chosen by GF03 
and U03 respectively), and Ecut needs to be reduced to a relatively 
narrow range of values around 110 keV in order to match the XRB 
up to 100 keV. 

The large or = 0.5 required above, however, implies that the 
number of sources with an intrinsic spectral slope harder than that 
of the XRB (F = 1.4) must be larger (15 per cent) as compared to 
the case of or = 0.3 (5 per cent), an increase by a factor of three. 
There is no obvious evidence for such a population. Many X-rays 
surveys with enough count statistics to be able to independently fit 
the AGN intrinsic power-law and cold absorption components have 
found a mean intrinsic p hoton-index <T>k, 1.9, with a disper- 
sion of ~ 0.2 — 0.3 (e.g., Nandra & Pounds 1994; Piconcelli et al. 
20031; IPerola et al J [2 004 ; Mateos et al. 2005b; Tozzi et al . 2006; 

independent of redshift. IPage et al.l j2006l) 



Mainieri et alj|200( 



have found a slightly larger dispersion or ~ 0.4 but with a slightly 
larger mean photon-index (<r>= 2.0) as well, resulting in a com- 



parable fraction of hard-slope sources (7 per cent) to the first case 
above. Luminous radi o-loud AGN are known to have harder slopes 
than radio-quiet ones teeeves & Turnej2000l) . but the vast majority 
of sources constituting the XRB are radio-quiet (e.g. iBarger et al.l 
.2001 ). 

Furthermore, the above requirement of a small -Ecut is also at 
odds with observations. While some A GN with -Ecut value s around 
100 keV (or lower) are known (e.g. iMolina et alj|200^ , studies 
with different satellites so far point to a much l arger mean cut- 
off energy, along w ith significant intrinsic spread jZdziarski et al.l 
ll995l ; lDadinall2007h . The on-going Swift/BAT survey does not see 
any deviation from an E~^ photon-index power law model in the 
100-200 keV band in the sum of AGN spectra observed, suggesting 
that any break in the spectrum should be at > 150 keV (J. Tueller, 
priv comm; see also iBeckmann et ai] |2007h . This is an important 
result that needs to be confirmed with the planned increase in the 
number of BAT channels analysed. 

It has been suspected that there may be a systematic calibra- 
tion uncertainty of the HEAO- A2 detector only, leading to an un- 
derestimate of the counts below ^ 10 keV, but leaving the XRB 
flux at and above its 30 keV peak unchanged, effectively 'steepen- 
ing' the XRB spectrum. This may then alleviate the need for signif- 
icant reflection, though the dispersion in F required will be close to 
the extreme observed values of or ~ 0. 4; moreover, a well-defi ned 
exponential cut-off is still required (see lRevnivtsev et"al]|2005l for 
recent discussion on the cross-calibration of the A2 detector). Very 
new results also suggest that the properties of AGN above 10 keV 
differ from those determined in lower ene rgy surveys (e.g. hard er 
intrinsic slope, higher XLF normalization; ISazonov et al ] |2007bO 
- it is crucial to see whether these hold up once more energy bands 
and signahnoise is accumulated; current constraints are also consis- 
tent with a canonical photon index and the presence of reflection. 

The above discussions have shown that not including reflection in 
SEDs that synthesize the XRB leads to a deficit of the modeled 
hard X-ray flux that cannot be easily compensated by modifying 
SED parameters or including additional populations. Our aim was 
not to unequivocally rule out exotic populations or constrain al- 
lowed ranges of variables by an exhaustive exploration of parame- 
ter space; we simply wish to point out the difficulty of reproducing 
the XRB in the absence of reflection (especially in the light of the 
recent Swift results mentioned above), while at the same time sat- 
isfying the well-established constraints below 10 keV. For the rest 
of this paper, we include a reflection component in (at least some) 
AGN spectra, either due to light bending alone or in combination 
with a PEXRAV component with -R = 1 (in order to compare with 
previous studies). 



3 X-RAY SPECTRA OF LIGHT-BENDING MODELS 
3.1 Unabsorbed spectra 

We first compare AGN spectra with and without the effects of light 
bendi ng included (Fig.[^. The 'PEXRAV' (Magdziarz & Zdziars^ 
Il995h SED is taken to be the unabsorbed template SED of GF03: a 
power-law with a photon-index (F) of 1.9 and a reflection compo- 
nent from a flat, infinite slab covering 27r steradians viewed at an 
inclination angle (i) of 60 degs. The only difference is the exponen- 
tial power-law cut-off (Ecut) which has been increased from 360 
to 500 keV, in order to directly compare with the results of U03. 
The spectra that incorporate light-bending are plotted for a range 
of source heights (S06 generate spectra at four discrete heights of 
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Figure 1. (Top) Model XRB spectrum without any reflection included (to- 
tal: thick solid black line; split according to obscuring co lumn: thin solid red 
line s) compared to the o bserved background data (blue; iGruberet all 19991 
an d iBarcons et alj200A scaled according to U()3). The model with the con- 
tribution of Compton-thick AGN increased by a factor of 4 is shown as 
the dashed lines. (Bottom) Including a gaussian dispersion (crp) of photon 
indices in SEDs without any reflection. Solid li ne: same as figur e above; 
dotted line: or = 0.2 and Ecut=200 keV (c.f. lGilli et ai]|2007h : dashed 
line: Increasing crp to 0.5 hardens the 1-10 keV slope to approximately that 
of the XRB, but overpredicts the hard X-ray flux; red dot-dashed line: i?cut 
needs to be decreased to fa 110 keV in order to match the XRB peak and 
subsequent fall-off. 



h = 2, 5, 10 and 20). For the present comparison, these spectra 
have been averaged over all viewing angles, and are normalized at 
1 keV. Identical intrinsic T and iScut values are assumed in all the 
spectra plotted in the figure, including the reflection-free power-law 
('PL'). 

As source height decreases, more photons are bent towards the 
central object, and also impinge on the accretion disk, resulting in 
a larger reflection fraction. For the lowest height plotted {h — 2), 
the peak monochromatic flux is higher than that of the PEXRAV 
SED by a factor of 26 (6), for a normalization at 1 keV (over 2~ 
10 keV). As h increases, the photon source is located in regions 
of shallower gravitational potential and the effect of light bending 
decreases until, for h = 20, the reflection fraction closely matches 
the PEXRAV 27r reflector model. 

The result of an increased reflection component is that the 
spectra become more 'peaky' in EFe{= vF^) units. This is il- 



Figure 2. Comparison of angle-averaged broad-band light-bending SEDs 
(source height as labelled) with a PEXRAV (i? = 1; = 60 degs; solid 
black) and simple cut-off power-law (PL; no reflection; lowest dotted black) 
SEDs. All spectra assume an incident power-law with F = 1.9 and an 
exponential cut-off i?cut = 500 keV, and have been normalized at 1 keV. 
The observed XRB fit from lGruber et alj h999ll . also similarly normalized, 
is shown as the blue long-dashed line. 



lustrated in Fig.[3] where the low {Ei) and high (iJu), or upper, en- 
ergies corresponding to the full-width-at-half-maximum (FWHIVQ) 
for each spectrum are plotted as the left-most and right-most locus 
of unfilled circles. The FWHM for the h — 2 and h = 20 spectra 
are ~ 150 and 270 keV respectively. In terms of the more familiar 
log-log plots for these spectra (c.f. Fig. [2J, Eu lies one dex in en- 
ergy above Ei for the h = 2 spectrum, and 2 dex above Ei for the 
h — 20 spectrum. 

Also plotted is the locus of peak energies as a function of 
source heights averaged over all angles (filled circles), as well as for 
the cases of low inclination angles (averaged over 0-50 degrees; red 
squares) and high inclination angles (averaged over 50-90 degrees; 
blue squares). Assuming strict correspondence with the orientation- 
based unification scheme, the latter two should correspond approxi- 
mately to viewing angles appropriate to Seyfert Is (hereafter Sy Is) 
and Sy 2s, respectively. Furthermore, at low source heights a large 
proportion of photon geodesies intersect inner portions of the accre- 
tion disk, with the result that reflected photons can gain energy via 
relativistic aberration and frame-dragging, especially at large incli- 
nation angles (see S06 for details). For the case with h — 2, the 
peak rest-frame energy can lie above 40 keV. We finally note that 
the fluorescent Fe Ka line self-consistently included in the code of 
S06 is highly smeared when relativistic distortions become impor- 
tant (Fig.[2j. 



3.2 Absorbed spectra 

The sources identified in recent 2-10 keV surveys, and responsible 
for the hard slope of the X-ray background, are mainly obscured 
AGN. In Fig. [4] we show the effect of obscuration on AGN incor- 
porating light-bending (the case with source hei ght ft = 2 is shown 
for illustration). We use the Monte Carlo code of IWilman & Fabiai] 



^ These FWHM values correspond directly to widths measured at half of 
the peak flux; i.e., they are not from gaussian approximations. 
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Absorbed spectre with h = 2 
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Figure 3. Compaiison of peak energies and energies coiresponding to the 
width at half-maximum of spectra with different source heights. The peak 
energies are plotted as filled circles while the lower half-maximum (Ei) 
and higher half-maximum (Eu) energies are plotted as the empty circles 
at each h. The peak energies of Sy Is (i.e. averaged over small inclination 
angles of 0-50 degrees) and Sy 2s (50-90 degrees only) are shown as the red 
and blue dashed lines respectively, though without the effects of absorption 
included for now. The equivalent values fo r the PEXRAV model are plotted 
as hexagons at the top; and those for the iGruber et al.l h999h functional 
form fit to the extragalactic XRB are marked as the three green dotted lines. 



Figure 4. Spectra producing by the light bending model (in this case, with 
the source height h = 2), absorbed by different amounts of absorbing 
columns of gas (logA^jj) labelled. 



l ll999l) , as implemented in GF03, to fully simulate the effects 
of Compton down-scattering (including the Klein-Nishina cross- 
section), which are especially significant at high column densities. 
Fig.|5]shows the widths and peak energies of these absorbed spectra 
(h = 2). As expected, large column densities result in further in- 
crease in the sharpness of the spectra, with the highest column den- 
sity of log7VH=24.75 resulting in a FWHM decrease by a factor of 
three, as compared to the unabsorbed spectrum. At high columns, 
the decrease in the peak energy due to Compton down-scattering is 
also evident. For the models considered in the following sections, 
spectra at intervals of AlogA'^H=0.5 ranging from logA'^H=20-22 
(Sy Is; intrinsic spectrum averaged over viewing angles of 0-50 
degs) and logA'^H=22-25 (Sy 2s; averaged over 50-90 degs) are 
created, for each height h. For logA'^H>25, photons over the entire 
X-ray range are severely depressed and we neglect the contribution 
of such sources to the XRB (e.g. iWilman & Fabiaa,1999,) . 



4 X-RAY BACKGROUND MODEL 

We now investigate the extent to which light bending spectra can 
contribute to the cosmological XRB. In order to do this, we re- 
quire a model of the distribution of AGN and their evolution. 
While recent studies with Swift and INTEGRAL have begun to 
probe large samples of AGN over the energy range corresponding 
to the peak of the XRB, these are still sensitive mainly to the lo- 
cal p opulation, and resolve at most a few pei* cent of the relevant 
XRB ( Markwardt et al.l 20051 ; Iseckmann et al.l2006l ; ISazonov et al.l 



l2007b ). On the other hand, the Compton-thin population of AGN 
is detectable below 10 keV, and has been extensively studied over 
the past few years, from which it is possible to extrapolate to higher 
energies with the SEDs of the previous sections. 

One recent comprehensive determination of the 2-10 keV X- 
ray luminosity function and A^h absorbing column density distri- 




E (keV) 



Figure 5. Peak energies and widths of the spectra from Fig.|4](flare height 
h = 2) absorbed by different columns of cold gas. See Fig. [5] for other 
details. 



bution is that of U03. A luminosity-dependent evolution in the 
nu mber density (LDDE) was inf erred by these auth ors and also 
bv lLa Franca et alj ( I2005L but see lBarger et alj|2005l) . These stud- 
ies also infer that the A^h distribution fraction (equivalently, the 
ratio of Type 2 : Type 1 AGN) probably evolv es with luminos 



ity an d/or redshift (see also Hopkins et al. 2007, Bal lantvne et al 



loot ; and for opposing vie ws: .Perola et al.. .2004,., Wang & Jians 



20061 iDwellv & PagellioO^ . Since our purpose is to investigate 
the perturbation on the XRB spectrum caused by inclusion of light 
bending in the spectra, we will simply use the LDDE model and 
A^'h function of U03 as our baseline for the distribution our AGN; 
similar relative modifications will apply in the case of other XLFs. 

Given the uncertainty in the population size of Compton-thick 
sources (U03 directly determined the number density of Compton- 
thin AGN and used, as a rough a pproximation f or Com pton-thick 
sources, the A^h distribution of iRisahti et all il999h for local, 
optically-selected Sy 2s in order to reproduce the total XRB spec- 
trum), we allow the fractional contribution of Compton-thick AGN 
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(logA'^H=24-25) to vary by a factor of ^ rcT ^ 2, where 
rcT = 1 represents an equal normalization between Compton- 
thick sources and Compton-thin ones with logA'n =23-24, as as- 
sumed by U03. Finally, in order to account for uncertainties in the 
XLF (at least 5 p er cent) as well as the shape and normaliz ation 
of the XRB itself jChurazov et"ai]|2006l : Iprontera et al.ll2007l) . we 
allow a final variable normalization (r) of between 0.95 and 1.05 
applied to the overall spectrum. 

Fits to the XRB are carried out over the 3-100 keV range by 
assuming that the HEAO error bars are normally-distributed, and 
computing a weighted sum of square deviations (an effective 'x^')- 
The difference in this fit parameter (Ax'^) is computed relative to 
the base U03 model (without inclusion of light bending). Given that 
observations from several different instruments are used, and since 
U03 do not actually fit their model to the observed XRB, such a 
approach may not be strictly valid. Due to the very distinct prop- 
erties of the light bending SEDs that we investigate, however, such 
a fit parameter will be a sensible measure of how well a model 
describes the data, at least relative to the base model. Our inten- 
tion is not to devise a new synthesis model of the XRB, but rather 
to explore the limits that can be placed on the prevalence of light 
bending SEDs within the existing constraints. We thus simply look 
for solutions that give a 'comparable or better' fit (Ax^ ^ 0), and 
consider all these to be acceptable. 



5 RESULTS 

Any realistic AGN population is likely to have a mixture of sources 
that readily show light bending, and those that do not. We can rea- 
sonably expect the heights of individual flares (on typical scales of 
a few Rg) to be distributed independently of the obscuring column 
density (on scales of at least thousands of Ro, if due to the torus). 
Furthermore, there is no apriori reason to assume any particular 
redshift-dependence in the reflection fraction of AGN, since X-ray 
spectral properties and BH spin distributions seem t o extend from 
their locally observed values out to high redshift (e.g. IVignali et al.l 
[2003; Volonteri et al. 2005). 

We thus assume that light bending is important in a certain 
fixed fraction (/) of all AGN, irrespective of obscuring column den- 
sity, luminosity and redshift. SEDs used to synthesize the XRB are 
generated as a simple linear combination of light bending templates 
and either i) PEXRAV SEDs (§|5TT); or ii) simple cut-off power-laws 
(i^ 15.21 ). The SEDs are normalized to have the same 2-10 keV un- 
obscured luminosities, combined in the ratio f ■ 1 — f (thereby 
conserving the number density of sources as given by the XLF) and 
obscured with various column densities of cold gas. For simplicity, 
we investigate several discrete values of / only (though covering 
the full range 0-1). Table [T] lists, for each height h and fraction 
/ investigated, the solution with the best-fit rescaling factors for 
Compton-thick sources and for the overall spectrum, as detailed 
below. 



5.1 PEXRAV + light bending scenario 

Firstly, in order to compare with previous XRB synthesis models, 
all base template spectra are taken to include a PEXRAV reflection 
component with R — 1, and a fraction / of these additionally 
include a light bending component. Fig. [6] shows the acceptable 
models (in terms of Ax^) as a function of percentage contribution 
(/ X 100) at various flare heights. 

The case where extreme light bending dominates in a large 




Percentage contribution (= fxlOO) 



Figure 6. Ax^ for the acceptable models incorporating light bending with 
different percentage contributions (/ X 100) of light bending. The change in 
is computed with respect to the base PEXRAV model with no light bend- 
ing included. Empty circles represent solutions in which the best-fit fraction 
of Compton-thick sources (rex) is small (< 0.5). In these cases, a sec- 
ondary fit was carried out by forcing 0.5 ^ rcT ^ 2 (boxed circles), which 
results in the apparent jaggedness of the othei'wise-smooth Ax^ curves. 



population of AGN is strongly rejected. For h = 2 (blue dotted line 
in Fig.|6ll, the resultant spectra are very 'peaky' and only up to ~ 2 
per cent of all AGN may include such pronounced light bending, 
if the XRB is not to be overproduced; the solution with /=0.05 
is also possible, but requires the near-total absence of Compton- 
thick sources (rcT=0.03; see Table [T) in order to compensate for 
the peakiness of the SEDs. In contrast, a significant population (up 
to /=0.1 with rcT>0.5, or /=0.3 with rcT=0) of all AGN could 
have a contribution from light bending with flares lying at ft = 5 
(black solid line). Finally, a fraction as large as 50 per cent of all 
AGN could include a significant contribution from light bending 
with h — 10 (red dashed line); these sources are consistent with 
reflection fractions (ii) of value 2 or larger (see S06). A fraction 
f=l is also allowed, but only with r-cT=0. Since the h — 20 SEDs 
match the base PEXRAV SEDs well (Fig.[2ll, we do not investigate 
h — 20 solutions further for now. 

The larger range of allowed fractions of h — 10 SEDs as 
compared to the range at lower h simply reflects the smaller SED 
peakiness. Several solutions (at all h) prefer low values of rcT (Ta- 
ble[Tll due to the fact that the peaky nature of Compton-thick AGN 
(either with, or without light bending) can be mimicked by light 
bending SEDs. Since evidence suggests that a sizable fraction of 
Compton-thick AGN is likely to exist, an additional fit was carried 
out for the solutions with rcT< 0.5 by fixing 0.5 ^ rcT ^5 2. For 
all bad fits, another solution with slightly looser constraints on the 
overall normalization r was searched for. In any case, all accept- 
able solutions have r values lying between 0.93-1 (Table[T). As an 
additional test of the quality of the fit, a simple power-law was fit 
to the total modeled XRB over 3-10 keV. All acceptable solutions 
over the 3-100 keV range match the base U03 model 3-10 keV 
slope (r = 1.35) to within | Ara-iol < 0.03. 

The overall fit statistic deviation of the acceptable solutions 
tends to decrease with increasing contribution of light bending. At 
large heights (say h = 10), this could be attributed simply to the 
presence of two reflection components, which may be associated 
with a disk-l-torus system (say). At lower heights, the stronger ef- 
fects of light bending responsible for the better fits (an increase in 
the reflection and peak energy of the overall modeled spectrum due 
to higher fraction of high energy photons present) are more appar- 
ent. This is better illustrated in the following section. 
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Table 1. Results of the fits for the models where a fraction / of sources include light bending, all with flares at height h. Ax^ is computed with respect to the 
U03 model without light bending. Typical 90 per cent errors (Ax^=2.71) on the Compton-thick (^ct) and overall (r) renormalization factors are 0.03 and 
0.01 respectively. An asterisk (*) indicates models where the minimum value of rex was fixed to be 0.5 during the fit. Fs—io is the fitted power-law index 
for the 3-10 keV modeled XRB. Solutions for the two scenarios discussed in §|5]are tabulated - PEXRAV + light bending: ^ 15.11 PL + light bending: 15.21 
At each h, the acceptable solutions (Ax^ ^ 0) with the maximum allowed values of/, and which also allow rcT>0-5, are listed in bold typeface, while all 
unacceptable solutions are in italics. Unacceptable solutions far outside the allowed range of / values are not listed. 



5.2 PL + light bending scenario 

We now investigate the possibihty that the only reflection compo- 
nent present in AGN spectra is that generated by light bending. 
The base spectra are thus simple F = 1.9, i5cut=500 keV power- 
laws not including any PEXRAV component, to which a certain frac- 
tion / of sources with light bending due to flares at height h, is 
added. Neglecting the PEXRAV component means that we also ne- 
glect sources of reflection other than the disk itself, e.g. the inner 
torus walls. This may be a simplification, but it is important to sep- 
arate the different reflection components and ask how much light 
bending alone can contribute, in order to clearly see its effect. 

Fig. |7] shows the acceptable solutions obtained by fitting the 
XRB over 3-100 keV. Only limited, well-separated ranges of 
fractional contributions of light bending are allowed for different 
heights. The more peaky the SED (i.e. the lower that h is), the 
less is the maximum allowed value of /, as was also found for 
the PEXRAV -^light bending scenario in the previous section. On the 
other hand, since reflection-free SEDs underestimate the XRB peak 
significantly (see §[2l(, very low values of / are excluded, in contrast 
to the previous scenario. Indeed, several solutions (TableO require 
r and/or rcT values 1 in order to compensate for the flux lost 



due to the absence of the PEXRAV reflection component. Solutions 
in the previous section all required normalization values of < 1. 

The middle plot in the same figure shows, as an example, the 
XRB residuals for one of the acceptable solutions (h=5; /=0.3). 
The peaky nature of the SEDs is able to reproduce the peak of the 
XRB in place of the PEXRAV component. Only comparatively mild 
renormalization factors (rcT=0.59±0.03 and r=1.00±0.005) are 
required in this case. Furthermore, the increased high energy (hard 
X-ray) flux of the light bending SEDs (Figs. H [3j can shift the 
XRB peak to energies closer to the observed peak (by about 10 
per cent in this case), and the residuals around 13 keV and 40 keV 
are reduced by a factor of ~ 25 per cent compared to the case 
with no light bending. The slight overestimate of the XRB between 
70-100 keV is not serious in terms of the fit statistic, and could 
ultimately by reduced by adopting a somewhat lower iJcut value 
(without affecting the spectrum at lower energies). 

The bottom plot shows the 2-10 keV source counts predicted 
by this model, showing that the U03 result is very well reproduced. 
Though the fraction of Compton-thick sources is reduced to 59 per 
cent as compared to the base U03 solution (leading to the slight un- 
derestimate of the U03 prediction [black dot-dot-dot-dashed line] 
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Figure 7. (Top) Ax^ for model fits with reflection-free power-law SEDs 
that incorporate different percentage contributions (/ X 100) of light bend- 
ing. See Fig.|6]for other details. (Middle) F ractional model diff erence spec- 
tram shown with respect to the XRB fit of lGraber et alj i 19991 , scaled ac- 
cording to U03), for one of the solutions (/i=5, /=0.3). (Bottom) 2-10 keV 
source counts predicted for this solution. 



at the faintest fluxes by 8 per cent), we note that such a reduction is 
not a characteristic feature of all these solutions (Table[Tll. 



6 CONSEQUENCES AND REFLECTIONS 

6.1 XLF based on intrinsic flare luminosities 

The U03 XLF involves the intrinsic 2-10 keV luminosity (L2-10) 
and inherently assumes that this can be computed by correc- 
tion of the observed source flux for obscuration and redshift ef- 
fects only. On the other hand, one of the main conclusions of 
iMiniutti & Fabiaij j2004h was that, under the light bending sce- 
nario, variations in the height of a flare of constant luminosity 
could cause the observed primary power-law emission to vary 
dramatically (while leaving the reflected component largely un- 
changed). For properly 'flux-calibrated' spectra determined by S06 
and shown in Fig. [S] (see also Figs. 6a and 15 of S06), the flux 
seen at infinity decreases with decreasing flare height (with a de- 
pendence on the inclination angle) due to multiple reflections, scat- 
tering, gravitational redshifting and photons falling into the event 
horizon. 

Light bending thus necessitates a distinction between the lumi- 
nosity inferred at infinity (Loo) without accounting for such losses. 



and the primary flare luminosity (Lflaro)- The modeling in the pre- 
vious sections has assumed the former case (L = Z/2-10 = Loo), 
since the U03 XLF (d"I>/dlogL) has been determined from the ob- 
served flux, corrected for A'^h and z. On the other hand, if one 
knows the fraction of intrinsic flux from a flare that is scattered 
out of the line-of-sight, it should be possible to recover the intrin- 
sic XLF based on Lflarc- Such an estimate is, of course, model- 
dependent but nevertheless useful in order to see the overall impact 
of light bending, since the AGN accretion power should be ulti- 
mately related to the intrinsic flare luminosity itself. 

For the fraction / of sources in which light bending is im- 
portant, let foo be the ratio between the 2-10 keV luminosity in- 
ferred at infinity to the intrinsic 2-10 keV flare luminosity (Loo ~ 
/ooLflarc)- Then, 



dlogLfla 



= (1-/) 



+ / 



dlogL • dlog(L/o> 



(1) 



$ is a function of L and z, and in the investigations of the pre- 
vious sections, we have assumed / to be constant, except for the 
dependence on flare height h. 

The best way to determine foo is to follow the energies of indi- 
vidual photons on their trajectories, including reflected geodesies, 
but this is beyond the scope of the present work. Instead, we com- 
pute foo approximately by using the flux calibrated SEDs of Fig. [8] 
and using the fact that the intrinsic flare spectrum will have the 
same shape as the primary cut-off power-law component in the light 
bending SEDs. Were light bending effects unimportant, an isotropic 
flare would have half of its primary flux reflected off the disk, and 
half observed directly at infinity; this is approximately true for the 
h = 20 spectrum (Fig. [8j. Thus, the fraction of the intrinsic flare 
luminosity that is observed beyond the sphere of influence of the 
BH's strong gravity can be determined by integrating over the spec- 
tra in Fig.[8]relative to the integrated 2-10 keV flare spectrum. We 
call this fraction /h//flarc, and list the measured values in Table|2] 

Another correction is required for the rest-frame flare flux lost 
down the BH event horizon (/lost), the cross-section for which 
increases with decreasing source height. This can be determined 
through simulations by computing the fraction of photon geodesies 
emanating from a height h that intercept the event horizon. Thus, 



Loo[h) — Lfla 



Jk_ 

-^flarc 



(1 - /lost) 



giving the required fraction (foo) for Eqn.[T]as 

Ih 



foo{h) 



/fla 



X (1 - /lost) 



(2) 



(3) 



The above correction fractions are listed in Table|2l 

Fig. m shows ci$/cilogLflarc (Eqn. [T} computed by using the 
maximum values of the light bending contribution (/; Table \T\ 
bold) for the 'good' (rcT>0.5) solutions of Figs.[6]and|7] in con- 
junction with the above values of foe- We choose 2 = 0.1 simply 
for illustration (this is the lowest redshift plotted in Fig. 1 1 of U03). 
The number density of sources (per cubic comoving Mpc) can ex- 
ceed the XLF without light bending by a factor of between ~ 5 and 
160 for the PEXRAV -Flight bending case, and between 7 and 770 for 
the PL-l-light bending case; the largest excess being at highest lu- 
minosities due to the steep slope of the XLF there (relevant for the 
second term on the right in Eqn.[Tll. Since / is taken to be indepen- 
dent of z, similar behaviour is expected at higher redshift, with the 
exact deviation from the base XLF depending on its slope. 

The radiative Bolometric luminosity (Lboi) of AGN is deter- 
mined by summing the observed SED over all wavelengths (with 
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the big blue bump [BBB] emission, usually reprocessed into the in- 
frared regime, dominating the output). Under the light bending sce- 
nario, the intrinsic 2-10 keV flare luminosity is reduced by a factor 
of /oo • Assuming that part of the BBB luminosity is similarly scat- 
tered out of the line-of-sight, LboI itself will need to be corrected 
upwards by f^^. The Eddington fraction would then also be cor- 
respondingly increased (for a constant black hole mass). On a cos- 
mological scale, using the results of Fig. |9] would imply that local 
Seyferts should have mean Eddington ratios larger by a maximum 
factor of ~ 1.5 — 6 (when averaged over all logL < 44.5), while 
the corresponding increase for higher luminosity quasars should be 
- 5 - 500. 

The above increase in mean Eddington ratios for Seyferts is 
comparatively modest, and has already been inferred for at least 
one source (Miniutti et al. 2007h) . The increase for quasars, how- 
ever, is very large, especially since quasars are already thought to 
be accreting at rates close to the Eddington limit. Of course, since 
we have used the maximum allowed values of / from Table[T]in the 
above computation, more moderate (intermediate) excesses of the 
XLF and of the Eddington ratios are more likely. Furthermore, the 
largest excesses are due to the lowest flare height of h = 2 alone. 
Considering the higher flares at h = 5 and 10, the maximum Ed- 
dington ratio increases are smaller (1.5-2.6 and 5-14 for Seyferts 
and quasars respectively). Such an aposteriori comparison could ul- 
timately help to constrain the universal distribution of flare heights 
themselves. 

On the other hand, the above comparison of the XLFs may 
be telling us that light bending effects ought to be less prominent 
for highly luminous sources. It is already known that reflection de- 
creases as a function of luminosity, but this effect has only been 
observed for the neutral component of reflection o ff matter dis- 
tant from the AGN itself ( Iwasawa & Taniguchi 1993). Whether a 
similar decrease of the disk reflection component in the innermost 
regions of the accretion flow is real or not will require not only 
exten ding the current surveys to a larger sample jOuainazzi et al.l 
l2006l) . but also prope r modeling of ionized re flection, as expected 
for luminous sources teeeves & Tumeill2000h . 

Finally, a likely alternative is that LboI could remain largely 
unchanged inspite of the increase in L^^^ . This is because the 
BBB emission is attributed to thermal emission in the accretion 
disk, and the fraction of this radiation bent due to strong gravity 
need not be the same as that from flares, especially if a stress-free 
inner boundary condition forces the disk emission to peak further 
out from the disk. If LboI >> ^flaro" remains largely un- 
changed, then it is the Bolometric correction factors used to scale 
from observations over finite bands (e.g. ~ ^Boi/iobs^" 

over 2-10 keV) that would have to be correspondingly larger, 
as compared to the case where light bending is not important 
and all the flare luminosity were observable at infinity (k^~^^ = 

6.1.1 Increase in BH mass density due to radiation lost down the 
event horizon 

Since sizable light bending fractions seem to be allowed from our 
results in previous sections, and anywhere between ~ 1 ~ 23 per 
cent of the emitted power in these sources falls directly down the 
BH (Table [2}, it is natural to ask if the resultant increase in BH 
mass is also significant. One cosmology-free method of estimat- 
ing the total BH mass density accreted over the history of the Uni- 
verse is by correcting observed radiation background intensities in 
any given energy band for the effects of absorption and scattering. 
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0.023 
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0.02 
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0.180 
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0.17 
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10 


0.363 


0.011 


0.36 


0.17% 


0.27% 


20 


0.524 


0.003 


0.52 







Table 2. Fraction of 2-10 keV flux from a flare at height h observed at 
infinity (Ihl Ins-rc) and lost down the event horizon (/lost)- Simulations 
of 0.5 million photons isotropically emitted from an on-axis point source 
at height h were used in conjunction with the flux calibrated spectra of 
Fig.[8]for this computation, foo can exceed 0.5 sUghtly because the disk is 
truncated in the simulations. AA/n, is listed as a percentage of the usual 
total mass buildup A/bh (without light bending), for radiative efficiency 
Tj = 0.1. The last two columns refer to the scenarios - ^: PEXRAV+light 
bending (§113); ^: PL+light bending (§|0- 




E (keV) 

Figure 8. Light bending spectra with different source heights, and proper 
relative normahzations due to flux scattered out of the line-of-sight in the 
strong gravity regime. These spectra have been averaged over all viewing 
angles. The top solid, blue line shows the intrinsic flare spectram assumed 
[though without correction for flux lost down the BH, which will push the 
flare spectrum higher by a further factor of [1— /lost (^i)] ^^], with an identi- 
cal power-law photon-index and cut-off as the other spectra. The rest-frame 
2-10 keV range is de-lineated by the vertical green dotted lines. 



The total rest-frame radiative energy density -Brad can then be in- 
ferred by assuming an appropriate Bolometric correction factor for 
that band, and knowing the average source emission redshift. If this 
represents a fraction r; of the total accretion energy, the mass den- 
sity will be AfBH = (1 — r;)/?7 x Era.d/(? ■ Such a procedure was 
proposed by S oltaij jl982l) . a nd has been impro ved upon since by 
[pabian & Iwasawal 1 1999) and E lvis et"ai] ( l2002h . among others. 

The mass gain due to light rays bending towards and falling 
into the BH (AA/ib) can be estimated in a similar maimer, except 
for the fact that the radiative efficiency (if) is not involved because 
the mass increase is directly due to infalling radiation only. As 
we saw in the previous section, the inferred luminosity at infin- 
ity needs to corrected by a factor /oo (in addition to cosmological 
fc-correction factors) in order to derive the intrinsic flare luminosity 
for a fraction / of sources in which light bending is important. A 
fraction /lost of this is lost down the BH (cf. second term on the 
right in Eqn. [3}. Since / is taken to be independent of L, AA/ib 
can be computed relative to the total mass density accreted without 
correction for light bending (Mbh) simply as 
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Figure 9. The 2-10 keV XLF from U03 (black solid line), i.e. without cor- 
rection for light bending, and with L2-10 = Lod- The other lines show 
the XLF based on the flare luminosity (L2-10 = ^flarc). as derived from 
Eqn.[T] and using the results of the maximum / values derived at each h 
for the i) PEXRAV + light bending scenario (Left; cf Fig. |6), and ii) PL + 
light bending scenario (Right; cf. Fig.|7). The redshift of 2 = 0.1 is simply 
chosen as the 'local' example (cf. Fig 1 1 of U03). 
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Assuming a typical value for the radiative efficiency of r; = 0.1 and 
using the maximum values of / for the good solutions in Table [T] 
we find AMib/MaH values of less than 1 per cent (see Table O, 
except for the lowest source height {h = 2) where the corrections 
and /lost are largest. The case with / = 0.1 (PL+light bend- 
ing) being the extreme with AMvo/Mbh = 14 per cent; for refer- 
ence, the base A/bh accreted in this case is 3.9 — 9.7 x lO'' 
Mpc""^ for a Bolometric : 2-10 keV luminosity correction factor of 
20 - 50. 

We note that only single reflections were treated in the compu- 
tation of /lost in Table|2j the tabulated values are thus a lower limit 
due to the reprocessed/reflected disk photons that are captured by 
the BH. The correction for this is, however, not expected to be very 
large. Also note that with ?7 > 0.1 (relevant for rotatin g BHs), the 
derive d extra mass loss will increase. For r] = 0.15 (e.g. lElvis et al.l 
I2OO2I) , AMib/MBH « 22 per cent for the above case with h = 2. 

To summarize, AMib lies below current observational lim- 
its on Mbh, and can be neglected on a cosmological scale, unless 
flares occur only in the deepest parts of the gravitational poten- 
tial wells (h — 2), where corrections for gravitational redshift are 
largest. As mentioned in the previous section, however, the above 
case with h — 2is likely to be the extreme scenario, and / can rea- 
sonably be expected to be less than the allowed maximum values 
(Table[TJ bold) that we have used, leading to a smaller AMib. 



6.2 Predictions for hard X-ray bands 

The reflection-dominated nature of the spectra incorporating light 
bending makes them good targets for future hard X-ray mis- 
sions capable of imaging above 10 keV. The proposed Japanese 
NeXT (New X-ray Telescope/Non-thermal Energy eXploration 
Telescope) satellite is the most sensitive of the such planned mis- 
sions with the shortest time-scale till launch (currently, 2012). 
While the full energy range of sensitivity will be from 0.5-300 keV, 



two hard X-ray supermirror telescopes will provide an effective 
area of between 500-1000 cm^ at 30 keV. This is the largest sensi- 
tivity for hard X-ray focusing optics to date (it is about 10 larger 
than that of Suzaku/HXT), and will prove crucial for studies of 
the non-thermal Universe. The soft 7-ray detector on NeXT is also 
likely to be the first X-ray polarimeter in space. See lTakahashi et al] 
( I2OO6I) for more details on NeXT. We note that the formation flight 
European mission Simbol-X (with a launch date around 2013) 
has very complementary proper ties to the NeXT mi ssion over the 
~ 0.5 — 80 keV energy range (.Ferrando et al.ll200^ . On a longer 
timescale, the European Space Agency's XEUS mission (with a 
nominal ene rgy range of 0.1-15 keV, and a goal of up to 40 keV 
at least; e.g. IParmar et al]|2006h will be crucial for characterizing 
such sources at high-redshift, in addition to enabling source identi- 
fication via the sensitive studies of (skewed) Fe lines. 

The NeXT hard X-ray effective area gives a continuum sensi- 
tivity of several xlO~^ photons s~^ keV^^ cm~^ in 100 ks, or a 
point-source flux limit of ~ 5.0, 1.7 and 0.6 x 10~^* erg s~^ cm~^ 
over 8-80, 20-50 and 10-30 keV respectively, assuming a power- 
law source with a photon-index of 2 (the exact sensitivity depends 
on the final telescope configuration parameters, which are to be fi- 
nalized in the near future). At these flux limits, 35, 39 and 68 per 
cent of the total XRB in the respective bands will be resolved, for 
the XRB normali zation adopted by U 03. If the original XRB nor- 
malization due to lGruber et al.l ( Il999h is closer to the correct value, 
the above percentages will be increased by a factor of 1.26. 

As an example, in Fig.[TO]we show (as the thick lines) the cu- 
mulative XRB expected to be resolved by NeXT in the 8-80 keV 
band, along with the corresponding cumulative number counts as a 
function of flux. For this illustration, we chose one of the good solu- 
tions with the maximum allowed / value (in this case, / — 0.3 for 
h = 5 for the 'PL-l-light bending' scenario; cf. the bottom plots of 
Fig.O in order to show the effect of light bending. The predictions 
with no light bending are shown as the thin black lines in Fig. 1101 
While the total source numbers match well between the two cases 
(as they should, since we are ultimately fitting to the XRB spectrum 
itself), the number density of Compton-thin sources (green dashed) 
is higher than the case without light bending (black dashed) espe- 
cially at bright fluxes, while the number of unabsorbed sources is 
lower (red dotted). This is mainly due to the increased peakiness of 
SEDs viewed at large inclination angles (Sy 2s), which means that 
these sources appear relatively brighter and can be seen to greater 
distances. On the other hand, the number density of Compton-thick 
AGN (blue dot-dashed) is reduced with respect to that in the base 
model (black dot-dashed) at faint fluxes due to the lower renormal- 
ization factor rcT=0.6; Table [T). Again, at bright fluxes, there is 
a relative increase in the number density of Compton-thick AGN, 
due to the fact that their peaky spectra make them brighter than any 
given flux threshold that applies to the base model. 

At fluxes of a few times 10^^^ erg s~^ cm~^, the numbers of 
Compton-thin (22 ^\ogNii< 24) and unabsorbed AGN (logiVH< 
22) match each other closely, in line with recent INTEGRAL and 
Swift results in slightly different hard X-ray bands. At the nominal 
NeXT flux limit of 5 x 10"" erg s"^ cm"^, Compton-thick AGN 
(with logiVn = 24-25) will constitute about 1 1 per cent of the total 
counts (13 per cent for no light bending). In comparison, Compton- 
thick AGN constitute only about 2 per cent of the total counts in the 
2-10 keV band. With an expected super-mirror angular resolution 
of ~ 30 arcsec, the above observation in the 8-80 keV band will 
not be confusion-limited at the NeXT flux limit. 
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Figure 10. (Top): The cumulative 8-80 keV XRB expected to be resolved 
by NeXT, as a function of flux. The width of the blue band near the top de- 
notes thejionnaliza^^ uncertainty in the observed XRB: between 1 and 
1.26 iUeda et al. 2003) X the background level inferred from the fit of 
iGruber e t al. ( 1999). (Bottom): Cumulative number counts as a function of 
flux in the 8-80 keV band. The solution with h = 5 and / = 0.3 from 
Fig. [T] is chosen for illustration. The thin black lines show the case with 
no light bending. The thin yellow vertical lines are drawn at representative 
100 ks and 1 Ms flux limits of 5 and 1.6 X 10~^* erg s~l cm~^. 



6.3 Reproducing the XRB without high absorption? 

Comparing the shapes of the SEDs in Fig. |2] with the shape of the 
XRB leads one to question whether it is not possible to reproduce 
the peak of the XRB only with unabsorbed, reflection-dominated 
sources. A decrease in flare height h increases the dominance of the 
reflection component relative to the observed power law component 
- an effect qualitatively similar to that produced by an increase in 
the obscuring column density. For sources in which the obscuration 
is measured in low signalmoise spectra, or inferred from count ra- 
tios in bands covering only limited energy ranges, it is possible that 
a hard observed spectral slope is due to light bending from flares 
at low height, rather than photoelectric absorption. So what is the 
XRB predicted under the radical assumption of using only light 
bended spectra without any absorption? 

We investigate this by using a crude parametrization as fol- 
lows. As before, we use the 2-10 keV XLF of U03 to describe 
the space distribution of AGN with intrinsic 2-10 keV luminosi- 



ties (Loo) inferred at infinity. But since the XLF is derived by 
absorption-correction of the observed hard spectral slope of AGN 
to some typical intrinsic power-law shape (F ~ 1.8 — 2), we should 
correct the light bending SEDs for some apparent absorption. As 
discussed in §[3] and ^ \6.U the unabsorbed template spectrum with 
h = 20 corresponds closely in shape as well as normalization to the 
unabsorbed PEXRAV case in which light bending is not important, 
and the isotropic flare flux can be inferred directly from observa- 
tions at infinity. Thus, we can use the h — 20 spectrum to compute 
the intrinsic luminosity that enters the XLF. Correction for 'appar- 
ent' absorption is then equivalent to normalizing the spectra at low 
heights h to the h — 20 spectrum in Fig. [8] 

The above correction is also approximately correct because 
the 2-10 keV fluxes of PE XRAV spectra absorbe d by columns of 
logiVH=21.5, 22.5 and 23.5 jWilman & Fabiarll9 99) are within 27 
per cent of the fluxes from unabsorbed light-bended spectra (Fig.jSj 
with h = 20, 10 and 5 respectively; while the flux of the Compton- 
thick logA'^H=24.5 PEXRAV spectrum is within a factor of 2 of the 
flux of the h — 2 unabsorbed spectrum. Thus, a simple association 
ofh = 20, 10, 5 and 2 with logA^H of 21-22, 22-23, 23-24 and 
24-25 will be correct, on average, to within the above uncertain- 
ties. For logA^H < 21, we can also use the h = 20 spectrum, which 
provides a good match over our energy range of interest 1 keV). 
With this association, it is then possible to also use the A'^h distribu- 
tion function of U03 to describe the the distribution of flare heights 
('/i distribution function') instead. 

Fig. [TT] shows the XRB spectrum that results from the above 
assumptions, i.e., using light-bended spectra with a distribution of 
flare heights h, with decreasing heights simulating the effect of in- 
creasing absorption due to cold gas, but without any extra obscu- 
ration included. Other than the excess below 3 keV, the agree- 
ment is excellent for the hard slope up to 10 keV, as well as for 
the XRB peak itself. In terms of the 3-100 keV fit parameters used 
in the previous sections, for this model, we have Ax^ = —122; 
r = 1.02±0.005; and rcT (now equivalent to rh=2)= 0.97±0.04. 
The excess below 3 keV can be removed by introduction of mild 
absorption (logA^n^ 22). 

The above result clearly shows the effect of unmodified light 
bending spectra alone. The good fit, without fine-tuning of the 
model parameters, is mainly due to the similarity between the 
'light-bended' and 'obscured PEXRAV' SEDs. We are not suggest- 
ing, however, that obscured AGN with logiVn >22 are unnecessary. 
Detailed comparison of the templates below 10 keV shows that 
photoelectric absorption can produce muc h steeper spectral ^radi- 
ents than even light bending at /i = 2 (cf. lWilman & Fabian" 19991 
and Fig. [8]l. Many recent surveys with high signalmoise spectra 
have indeed established the presence of significant obscuration in 
AGN environments, and we have used their results in the previous 
sections. 

On the other hand, the fact that the light bended spectra alone 
provide a good average description of the hard XRB slope above 
3 keV has an interesting consequence. It may provide another ex- 
planation for some AGN in which optical and X-ray classifications 
have been found to be mis-matched. In particular, recent surveys 
have established the presence of a population with apparently high 
absorption in X-rays, but with broad lines in the optical, implying 
no reddening due to d ust (^aiolino et al. 200ll ; IPerola et al.ll2004l; 
[Silverman et al.||2005l;rPage et al.ll200a). These constitute about 10 
per cent of the X-ray AGN population jGarcet et al .120071) ; the mis- 
match is attributed variously to the presence of dust-free gas within 
the sublimation radius, to large (optically-thin) dust grains, or to a 
dense ionized wind. Alternatively, in the light bending scenario, if 
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Figure 11. XRB spectrum predicted by the model with unabsorbed light 
bending templates and a distribution of flare heights. Decreasing h is ap- 
proximately associated with increasing logAfjj (§[63]and Fig. [8). There is 
excellent agreement with the observed XRB over ~ 3— 100 keV; the excess 
at smaller energies can be mostly removed by including mild absorption of 
logA^H < 22, as shown by the very thin black hne showing the total XRB. 



flare heights are predominantly low, a reflection-dominated (hard) 
spectrum will naturally result without the need to invoke obscura- 
tion; light bending will not affect emission from the optical broad- 
line region on scales of thousands of Gravitational radii. Studies 
find that optical-broad -line X-ray-absorbed AGN mainly have low 
column densities (e.g., iMateos et alj20053) . But, given current un- 
certainties in the selection effects of hard spectrum (and Compton- 
thick) AGN, the possibility that some of the fainter X-ray flat- 
spectrum sources (indicative of high obscuration) may have broad 
optical lines is not ruled out. How many of these mis-matched 
sources can instead be explained by light bending remains to be 
seen with the on-going follow-up of Swift and INTEGRAL de- 
tected AGN; the good match of the hard XRB slope that we find 
certainly allows for a significant proportion to be explained as such. 

6.4 The effect of ionization on reflection 

Since the ionizing flux of radiation at the inner edge of the accre- 
tion disk (rin, lying at the last stable orbit) may be intense enough 
to cause the surface layers of the disk to be ioni zed, it is important 
to comment on the effects of ionization (e. g., IPone et si 1 l 19921 ; 
IRoss & Fabiarilll993l ; iNavakshin et al.ll200(]|) . which has been ig- 
nored in our simulations. 

The latest ionized reflection model codes now include 
all relevant atomic transitions and ionization states (REFLION; 
lRoss&Fabiarill2005l) . Though the resultant spectra can be com- 
plex due to the host of emission and edge features produced, the 
major effect of ionization on the hard X-ray band of interest cover- 
ing the peak of the XRB is a change in the width of the Compton 
hump, with ionization parameters (^) larger than ~ 100 resulting 
in broader SEDs. For the canonical AGN population that is thought 
to contribute most of the XRB (Lx ~ W"* erg/s; A/bh ~ 10* 
M0), ^ can be estimated to have values of < 100 — 200 for flare 
heights (h) of 2-10 Rq- Here, we have defined ^ according to 
IRoss & Fabiani with disk density nu fixed at 10^^ cm~^ at a radial 
distance of TRg ■ This distance is where the emissivity of a standard 
(Schwarzschild) disk is expected to peak in the case of a stress- 
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Figure 1 2. Effect of Fe abund ance on ionized reflection spectra using the 
REFLION teoss & Fabian ''2005) code. The Ionization parameter g = 1000 
and photon index F = 1.9 were fixed. 



free inner boundary condition. The ionizing radiation flux will in- 
crease inwards, resulting in ^ values <^ 2000 at rin = 1.247?g for a 
constant density disk. Accounting for realistic disk density profiles 
will, however, ensure values of ^ smaller than this. It is also worth 
noting that in a recent Siizakii observation of MCG-6-30-15 which 
showed a steep inner emissivi ty profile, only a low ^ < 100 was 
inferred jMiniutti et"al]|2007ij) . 

In any case, for the above range of ^, abundance effects are 
likely to be more important than ionization. Fig. [T2] shows SEDs 
produced by REFLION for the relatively large, fixed value of ^ = 
1000, and with varying Fe abundances. Iron absorption over the en- 
ergy range of ~ 1—40 keV controls the steepness of the low-energy 
part of the Compton hump, as well as the peak hump energy. Higher 
abundances (as also inferred for MCG-6-30-15, for instance) result 
in peakier SEDs. 

These effects will be directly applicable to the light bending 
SEDs as well, with the Fe line of Fig. [12] being highly smeared 
in a Kerr metric. Ionized reflec tion is likely to be m ost important 
in 'Regime F identified by Miniutti & Fabiani ( l2004h when source 
heights are low and ^ is maximized. Searches with future missions 
for correlations between ionization effects (e.g. the ionized fluo- 
rescence Fe line and edge expected at ~ 7 keV, rather than 6.4 
keV from neutral matter) and light bending effects (e.g. the gravi- 
tational redshift of the line, or the strength of the Compton hump) 
may prove fruitful in th is regard. The blurring of sharp features due 
to relativistic rotation jCrummv et al.ll2006r) will, however, make 
this task non-trivial. 



7 DISCUSSION 

The hard X-ray background spectrum emphasizes energies of 30- 
40 keV, and seems to naturally suggest the presence of a reflection 
hump in sources that constitute the background. Indeed, in §|2l we 
discussed the difficulty in trying to reproduce the background spec- 
trum without any reflection, while simultaneously satisfying the 2- 
10 keV observational constraints. 

XRB synthesis models all incorporate some le vel of reflection, 
either by including a fraction attributed to the torus jGhisellini et al.l 
1 1994b . or simply assuming an infinite plane parallel optically-thick 
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slab reflector; in any case, an average value of the reflection frac- 
tion _R ~ 1 is able to fit the spectrum to within 10 per cent over 10- 
50 keV. On the other hand, there is recent evidence that, in at least 
some sources in the Universe, the reflection component is explained 
only by properly modeling photon trajectories that are bent in the 
strong gravity regime of the immediate vicinity of the BH. The net 
SEDs of these light bending models are predicted to be reflection- 
dominated and peaky within distances of a few gravitational radii 
from the BH (§ |3). Since the XRB is the integrated emission of 
AGN over redshift, it should be possible to use it in order to place 
constraints on such reflection-dominated sources. This is what we 
have tried to do in this paper. Our aim has not been to present a 
new XRB synthesis model, but rather to use the well-determined 
constraints derived from extensive 2-10 keV surveys, and to ex- 
trapolate to higher energies by using realistic SEDs in which light 
bending produces a char acteristically strong reflection component 
jSuebsuwong et al.ll2 006^. 

Assuming that light bending is important in a certain frac- 
tion (/) of sources, we derive broad constraints on / by fitting 
the hard X-ray background over 3-100 keV. / can be thought 
of as representing the distribution of sources with high effective 
reflection fractions {R > 1). Including both PEXRAV and light 
bending reflection components, we find values / ^ 0.02, 0.1 and 
0.5 (decreasing with decreasing source height h [or equivalently, 
with increasing spectrum peakiness]; i) l5.1t . Slightly larger values 
of / can be allowed if the peaky contribution of Compton-thick 
sources is neglected, though it is likely that a sizable fraction of 
Compton-thick AGN exists in the Universe. If light bending is the 
only source of reflection (the 'PL-l-light bending' scenario; § |5.2| l, 
then specific ranges of / « 0.05(/i = 2), 0.1 — 0.3(/!. = 5) and 
0.3 — 0.8(/i = 10) are required, due to the fact that reflection- 
free SEDs underpredict the XRB peak. Fi g . 1 1 3 1 summarizes the al- 
lowed solutions in the grid of models investigated (including now, 
for completeness, the solutions with h = 20, where SEDs are very 
similar to the PEXRAV templates). 

We derived the XLF based on the intrinsic flare luminosities 
(d$/dlogLflarc)> assuming the maximum allowed ranges of / al- 
lowed (§ I6.U . These XLFs were found to be higher than the case 
without light bending, by average factors of up to ~ 2 — 4 for 
Seyfert-like luminosities, and for flare heights h = 5 and 10. Av- 
erage Eddington ratios of local Seyferts may also increase corre- 
spondingly; the increase for sources with quasar-like luminosities 
is larger by at least a factor of two. Implications and caveats of 
this were discussed. The BH mass increase due to flux lost di- 
rectly down the BH was found be small, relative to the total mass 
computed without light bending corrections, except when h — 2 
(ii l6.1.1t . Realistically, / will probably be lower than the maximum 
values that we derive. Moreover, we have investigated light bend- 
ing from flares at only four discrete source heights, and assumed 
that only a single h value is dominant in all sources. The more 
likely, complex, case is that a distribution of all source heights will 
be present. Finally, we have only studied spectra that result from 
on-axis flares; for flares that lie off-axis, the high-energy flux can 
be further increased by doppler-boosting (S06). 

Since the Monte Carlo simulations of S06 assumed a spin- 
ning BH, how relevant are these SEDs for the XRB? Growing evi- 
dence ind icates that most super-massive BHs are probably rapidly 
spinning I Wang et al.| l206^; (eIvIs et alj|2002l ; iFabian et alj|2002bl ; 
IVolonteri et al. 2003), which means that the Kerr metric is probably 
a more relevant description of the space-time around BHs than the 
Schwarzschild one. This, in turn, allows for the possibility that the 
accretion disk extends down to the last stable orbit within a few Rq 



of the BH itself. Any photons reflected off the inner parts of such 
disks will be strongly affected (bent) by the gravitational potential. 
Thus, light bending should be considered seriously in a cosmologi- 
cal context. Our results show that the XRB allows for the existence 
of a significant fraction of sources with pronounced light bending. 

If the sources of X-ray photons are related to inner parts of 
this disk (either a reconnection flare due to magnetic lines that are 
anchored in the disk; or an electron corona evaporating off the inner 
portions off the disk), then these photon sources may lie close to 
the disk itself, at very low heights. Indeed, the highly-reflection- 
dominated sources in which light bending has been identified so 
far are consistent with flares lying within just a few Rq of an fast 
spinning BH (e.g. NGC 4051; S06), similar to the low h values that 
we have investigated. 

The increased hard X-ray flux of light bending SEDs viewed 
at large inclinations shifts their peaks to higher energies, and al- 
lows for slightly better fits to the XRB. Yet, this still leave some 
residuals in the XRB fits (over-predicting the HEAO spectrum at 
~ 15 keV by about 5 per cent, and similarly underpredicting the 
observations at ~ 40 keV in the best cases. This is a general prob- 
lem of all XRB synthesis models, since the rest-frame peak of the 
PEXRAV reflection component lies close to 30 keV. The rest-frame 
SED peak shifts to ~ 40 keV (Fig. [3] |5]l are not sufficient for 
an AGN population with an average redshift of z ~ 1, though 
the above residuals are smaller by ~ 25 per cent with respect to 
the case of no light bending. In this respect, super-solar metallic- 
ities could open up e xtra parameter spa ce by shifting the peak to 
still higher energies jWilman & FabianllT999,) . Such metallicities 
have been observed in th e spectra of many AGN (t he evidence for 
quasars is clear-cut re.g.. lHamann & Ferlandll999l l, while at least 
some examples exist for ca se of Seyferts as well [e.g., MCG- 6- 
30-15: lMiniutti et alj|2007ij and Mrk 1040: lRevnolds et al.lll995ll ). 
suggesting that a combination of high metallicity with li ght bending 
may not be unrealisti c. A population of ADAF sources ( ICaoll2007l ; 
ldiMatteoetal.ll999l) may also provide some contribution. Interest- 
ingly, we note that the XRB fits with ADAF sources of Cao (2003) 
also involved spinning BHs, with the best fits having a = 0.9. 

Another perspective of considering light bending is related to 
variability and state transitions. Transi tions in spectral shape s have 
been observed in several AGN (e.g. ICuainazzi et aLluOOS ^ and 
associated either with large changes of the obscuring gas column 
density, or with a transient dimming of the AGN activity. While 
timescales for such transitions are usually on the order of years, 
clear spectral variability spanning just a few days has been ob- 
served in at least one source (NGC 1365; iRisaliti et alj|2005l) . In 
the light bending scenario, a lowering of the height of the main 
emission region associated with flares of a given luminosity above 
the event horizon would also produce a transition to a dim, hard 
state, with this state being reflection-dominated. The flux of the re- 
flection component itself should remain roughly constant (as also 
observed in the case of NGC 1365). In such a scenario, our deter- 
mination of / is equivalent to determining the duty cycle of AGN in 
the reflection-dominated state. Again, our results allow for a sub- 
stantial fraction of sources to be in such a low, hard state at the 
present epoch. Strong light bending could then be the explanation 
for mis-matched optical and X-ray classifications in some AGN 
(§|63}. In such a scenario, a hard X-ray source spectrum would be 
due to the dim, reflection-dominated state, rather than obscuration. 

It is worth noting that our results do not include any sources 
that lie below the 2-10 keV flux limit of the deepest surveys con- 
sidered by U03. In this sense, our results may actually represent a 
lower limit to /. Since strongly light bended, intrinsically luminous 
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Ilowed contribution of sources with light bending vs source height 
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Figure 13. Allowed combinations of (/, h) investigated in § |5] At each 
height h, the lower, blue bar denotes the 'PEXRAV+light bending' sce- 
nario (Fig.|6) and the red, upper bar is for the ' PL+light bending' scenario 
(Fig. I?)- Acceptable solutions that require a small Compton-thick fraction 
(rcT<0-5) are outlined by the empty bars, the rest being denoted by the 
filled bars. Though bars are drawn continuously for clarity, only the dis- 
crete / values of the previous figures were investigated. 



sources will appear faint at infinity, they could contribute a com- 
paratively large proportion of the local mass accretion onto BHs 
without much contribution to the XRB intensity. Thus, correction 
for light bending in determination of the total local BH mass den- 
sity could also be used to constrain the prevalence of light bending, 
but since our approach has been to fit to the XRB its elf, we cannot 

[ lace useful limits via this method. Recent work bv lWorslev et al.l 
2005h has shown that the fraction of the XRB resolved in ultra- 
deep ,and XMM-Newton surveys decreases as energies approach- 
ing 10 keV are considered. They suggest that there is room for a 
population with an overall peaky spectrum at faint fluxes, and find 
that the missing fraction is consistent with a population of highly- 
obscured sources with logA'^n ~ 23 — 24.3 at 2 ~ 0.5 — 1.5. It will 
be interesting to see whether the peaky spectra incorporating light 
bending could also make up the missing population, though these 
would also probably have to be obscured if their redshifted flux be- 
tween ~ 5 — 10 keV is to lie below the threshold of the recent 
~ Ms surveys. Future hard X-ray imaging missions such as NeXT 
(ii l6.2t are in a good position to constrain the nature of these sources 
by sear ching for hard X-ray s from the X-ray undetected GOODS 
AGN of lWorslev et al.l ( l2006h . Precise large sky area measurements 
of the background near its peak will also prove crucial in order to 
constrain the contribution of light bending to the hard XRB. 

F inally, we note that on-going surveys with IN TEGRAL and 
Swift teeckmann et al.ll2006l:lMarkwardt et al.ll2005b should begin 
to detect sources in which light bending dominates, if these sources 
occur at low redshift and low-luminosity. Curre ntly, the fraction 
of Co mpton-thick AGN (e.g., in the catalogue of iBeckmann et all 
l2006t) is similar to that predicted in Fig. [TO] with a small differ- 
ence due to the bandpass assumed. Whether or not the hard X-ray 
spectra of detected sources can be distinguished from transmission- 
dominated spectra without light bending remains to be seen, once 
increased signahnoise is accumulated. Detailed study of the asso- 
ciated fluorescence Fe line profiles will help in this regard. 
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